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Abstract 



A catalogue of binary galaxies with 621 pairs has been determined by applying a surface 
density enhancement procedure to The Surface Photometry Catalogue of the ESO-Uppsala 
Galaxies. The method does not require any redshift information. An additional restriction, 
based on objective criteria that take into account the completeness of the source catalogue, led 
to a sample of 189 isolated pairs which are listed. We have obtained the optical luminosity 
function of binary galaxies in the catalogue, from which we estimate that the luminosity 
density of galaxies in binaries is ~4% of that found for field galaxies. The general properties 
of our sample are similar to those from CPG and CMC 

Key words: Galaxies: Binary — Galaxies: Clustering — Galaxies: Catalogue — Galaxies: 
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1 Introduction 



The clustering properties of the galaxy distribution in the Universe is one of the main 
topics studied by the observational cosmologists in the recent years. Clustering scale varies 
from "totally" isolated galaxies to clusters of clusters of galaxies. This interval is usually 
binned in different classes of multiplicity as a matter of simplicity, not necessarily implying 
any physical difference. In such hierarchy, pairs of galaxies represent the first step after the 
isolated galaxies with no clear cut between them. Their difference is only a matter of what 
kind of isolation criterion is used. However, even without a totally objective way to define an 
isolated pair of galaxies they can still be useful in shaping our understanding of the physics 
of galaxy interactions and formation mechanisms. 

The historical step in the way of understanding how interacting galaxies evolve and 
have their internal and orbital structure changed, came with the Toomre and Toomre's (1972) 
investigation on the M51 and Antennae systems. Modelling of interacting binary systems has 
now been done by many authors (e.g., Barnes 1988, Borne 1988, Balcells, Borne and Hoessel 
1989, Madejsky 1992, Keel and van Soest 1992). Close and wide binaries have been useful for 
the study of many aspects of galaxy morphology and its connection to galaxy formation and 
the Hubble sequence (Sulentic 1992, and references therein). Such systems are also useful for 
statistical analysis aiming the determination of galaxy mass (Karachentsev 1972, 1985, 1987, 
Turner 1976, Schweizer 1987, Soares 1989), and the study of the secular evolution of orbital 
parameters (e.g., Verner and Chernin 1987). 

Given the relatively large number of available binary lists, one might question whether 
the determination of yet another double galaxy list is really necessary. In fact, there exists 
a wealthy number of galaxy pair objective catalogues, i.e., catalogues determined with au- 
tomated and well-defined procedures. However, most of them covers the northern sky (e.g., 
Karachentsev 1972, CPG - Catalogue of Isolated Pairs of Galaxies), while the southern sky 
lacks of such an extensive survey. The available southern lists are either incomplete or based 
on subjective selection criteria (Arp and Madore 1985, Schweizer 1987, Zhenlong et et al. 
1989, Rampazzo and Sulentic 1992). 
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In section 2 we make a detailed description of the selection procedures used to build our 
binary sample. Section 3 is devoted to the determination and discussion of the luminosity 
function of galaxies in our binary catalogue. Finally, section 4 contains a brief discussion of 
the general characteristics of the sample. 

2 Sample Selection 

The selection of the present binary galaxy sample was based on the galaxy-galaxy 
projected separations and apparent magnitudes. The source list used for this purpose was The 
Surface Photometry Catalogue of the ESO-Uppsala Galaxies (hereafter ESO-LV, Lauberts and 
Valentijn 1989). The ESO-LV has 15467 entries, most of them with total magnitudes in B 
passband. The catalog is probably complete until 14.5 mag, as will be discussed in the next 
section and covers a solid angle of 2.01 sterradians. It is important to keep in mind that such 
incompleteness may introduce a bias in our binary sample, i.e., we may lose pairs for which 
the secondary galaxy is much fainter than the limiting magnitude. 

The likelihood of physical association for the candidate pairs was estimated following the 
prescriptions used in the construction of the Catalog of Multiple Galaxies (CMG, van Albada 
and Soares 1993, and Soares 1989; see also van Moorsel 1982, van Moorsel 1987, Oosterloo 
1988). A summary of the main points of such method is presented in the Appendix. For a 
more detailed discussion the reader is refered to Soares (1989). 

In what follows we show how the CMG method was applied to the ESO-LV. The 
method consists of investigating the frequency distribution of projected separations taking 
into account the apparent magnitude of each galaxy in the pair. This is done in the following 
way. Each galaxy is classified according to its apparent magnitude, and put in 0.5-mag 
bins, whose centers range from Bt,i = 12.0 up to Bt,g = 16.5. They are numbered by 
J = 1,2, ...,10. Thus, the interval J = n has all galaxies with magnitudes from Bt, u — 0.25 to 
Bt,u + 0.25. Now, a pair of galaxies is characterized by a given A J, defined as the difference 
between J(secondary) — J(primary). The distributions of separations for every AJ, from 
AJ = to AJ = 9, were investigated individually . We only considered pairs for which 
the primary galaxy was brighter than 14.5, minimizing the effect of incompleteness of the 
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source catalog in our final binary sample. However, the secondary was selected regardless its 
magnitude. Galaxies with no magnitude listed by the ESO-LV were included in the last bin 
(J = 10) and galaxies brighter than the 11.75 have been included in the first bin (J = 1) in 
order to prevent the undersampling at the brighter end. 

The distribution of first-neighbour projected separations is supposed to be Poissonian, 
in the absence of physical association. Any deviation, in excess over that, reflects small-scale 
gravitational clustering. The likelihood of physical association of a given pair is defined from 
the amount of deviation of the actual separation distribution over the random (Poissonian) 
distribution. Let p be such a probability, given by 



where x is the normalized projected separation (see below), pi{x) is the theoretical Poissonian 
distribution of x (fitted to large values of x). The observed values of x are obtained by dividing 
#i, the first-neighbour separation, by the Poissonian expected mean value of #i, which in turn 
depends on the surface density of galaxies in the vicinity of the pair (n). Such a normalization 
of Q\ , implies that the frequency distribution of x is independent of the local surface density 
of galaxies, making it possible to treat galaxies from different clustering environments into 
one distribution. Following the CMG recipe, we have used x as given by eq. (A5), 



in order to test for physical association. The actual distributions p {x) were obtained for pairs 
belonging to A J = 0,1,..., 5. We have adopted also a crude grouping of the morphological 
types in 3 broad classes, namely, early-type galaxies (—5 < T < 0), spirals (1 < T < 6) and 
late spirals and irregulars (T > 6). The excess of p {x) over the Poissonian distribution was 
fitted to a linear approximation for each morphological type of the primary. The solid lines 
shown in Figure 1 are given by 
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(2) 



5 



We can see from Figure 1 that the behaviour of p(x } T) is independent of the magnitude 
bin of the secondary, contrary to what was found for the CMG, where it is noted a substantial 
correlation of p(x } T) with A J (Soares 1989). Moreover, the relation is steeper for early type 
galaxies, since these objects are more clustered than the late type ones; such a behaviour was 
also noticed in the determination of the CMG. 

This parameter, p(x } T), was used as our primary indicator of association. If for a given 
bin of magnitude, p(x } T) is larger than a critical value (p c = 0.5) then the pair is selected 
as a probable candidate to form a binary system. In order to eliminate binaries in dense 
groups and clusters we have adopted a local density cutoff of 5 galaxies per square degree. 
This density measurement was obtained from the distance to the tenth nearest neighbour 
(see Appendix, or Soares 1989, for a complete discussion). Imposing all these constraints we 
have obtained a crude sample containing 621 binaries. Further examination has shown that 
70 binary pairs from this first list were actually members of multiple systems. This inspection 
was done looking for different pairs sharing the same galaxy. Excluding these obviously false 
pairs we got a purged list of 551 pairs. 

An important effect jeopardizing the binary selection process is related to the incom- 
pleteness of the ESO-LV catalog. A significant fraction of binaries in our sample has the 
secondary galaxy fainter than Bt = 15.0. For these pairs the separation tends to be larger 
than pairs with brighter secondary galaxies. We associate this effect to the incompleteness 
of the basic catalog resulting in a surface density of faint objects much less than what would 
be expected from a complete catalog in the whole magnitude range. The following procedure 
was used to minimize such effect. For a uniform distribution of galaxies the number of objects 
with magnitudes between m and m + Am can be expressed as AN oc 10°- 6m , and the mean 
separation between a primary and the nearest object is (8) = 1/2-^/n, where n is the surface 
density for galaxies with magnitudes in the previously defined interval. Expressing the mean 
separation in terms of the magnitude of the secondary galaxy we have (6) oc 10 _0 ' 3m , from 
which follows 
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log(0) = log(0 o ) - 0.3(m - m Q ). 



(4) 



This relation predicts the relation between 6 and m for a random sample. When physical 
association is present then the separation tends to be smaller than that predicted by eq. (4). 
Therefore, in a diagram 6 versus m binary sistems tend to locate below the relation expressed 
by eq. (4). 

Figure 2 shows the separation of the binaries in our sample, in arcsec, as a function 
of the magnitude of the secondary galaxy. Considering only the secondary galaxies in the 
interval 12.5 < Bt < m , where m — 14.5 is the limiting magnitude of the ESO-LV catalog 
we fitted the relation 

log(# max ) = 3.254 - 0.3(m - m Q ) 

The previous fitting was done taking for each bin of 0.5 mag the three largest separa- 
tions, and determining the upper boundary to the points. With such procedure we expect 
to have a fair representation of the separation of galaxies distributed uniformly as a function 
of the magnitude. The upper limit of the separation of galaxies in a uniform field represents 
the maximum separation binaries can reach on average, justifying why we substituted 6 by 
6 ma x hi the last equation. We have extrapolated the fitting function to m > m in order 
to eliminate those binary pairs for which m sec > m and \og(6) > log(# max ). Applying such 
restriction we were able to eliminate 106 pairs out of the 551, ending up with 445 binaries. 

We have then inspected a field of 2° X 2° around each binary in the sample, using the 
ESO-LV catalog to construct a map of the region. All pairs with probable companions inside 
a circle of 1 Mpc around the center of the pair were eliminated. We have used H = 50 
km/s/Mpc throughout this paper. Particular attention was paid to those cases where loose 
groups were easily identified in the outskirts of the 1 Mpc circle. Avoiding all these spurious 
cases we have defined our master list of binary sistems with 189 pairs, which are presented in 
Table 1. In each entry we list the NGC/IC number , or the ESO identification code, equatorial 
position for 1950 (RA and DEC), magnitude, numerical morphological type and velocity, for 
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the two components of the pair. Also listed is the separation, in arcsec, and the probability 
parameter p given by the equation 3. The last column lists the categorization of the vicinity of 
the pair based on the visual inpection of the maps. For those pairs with at least one available 
radial velocity we made a circle of 1 Mpc radius centered on the candidate pair. Two main 
categories were defined: (a) Probable pairs (P) - pairs with no companions inside the 1 Mpc 
circle, but with some suspicious companion in the neighbourhood. By suspicious we mean 
a galaxy with magnitude comparable to that of the pair's; (b) Isolated pairs (IP) - binary 
systems with no companions inside the IMpc circle, and also no clear companions in the 
neighbourhood. A probable pair with separation comparable to their component diameters 
was classified as CLP (Close Pair). Similarly for the isolated pairs, CLIP (Close Isolated 
Pair). We noticed that those pairs classified as isolated (IP+CLIP) are less frequent than the 
probable pairs (P+CLP). 

The presence of companions in the outskirts of the IMpc circle was coded as PCI 
(Probable Companion Outside 1 Mpc). Sometimes it is possible to see Faint Companions 
(FC) close to the primary, or secondary member. Finally, the presence of a nearby group 
outside the 1 Mpc region was coded as (NG). 

It is interesting to note that even discarding pairs in dense regions the number of pairs 
close to or superposed on larger structures such as loose groups represent 46% of the final 
sample of selected binaries. This percentage is fully consistent to what it is expected from 
simulations (see Mamon 1993). A cross check with ZCAT has shown that among the binaries 
with available redshift for both galaxies 20% are false pairs (AV > lOOOkm/s). This number 
can be regarded as an estimate of the contamination by spurious pairs in our sample. 

In Figure 3 we present the results concerning the completeness of our final sample based 
on the (V/V max ) test (Schmidt, 1968). Figure 3a shows the result for a sample of "field" 
galaxies extracted from the ESO-LV catalog, all of them with bn > 30°. By "field" we 
mean a sample of objects with known radial velocity and located in regions where the mean 
local projected density of galaxies is less than 0.7 galaxies per square degree, defining a field 
sample that is approximately 10% of the whole ESO-LV catalog. In Figure 3a filled and open 
symbols represent the field sample with magnitudes corrected and not corrected for Galactic 
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extinction, respectively. We can observe that in this case most of the "field" sample have 
(V/V max ) < 0.5. Applying the Galactic extinction correction of 0.12 cossec | b | most of 
the "field" sample have (V/V max ) > 0.5 and we can clearly see the gradual incompleteness 
starting at Bt — 14.5. Application of the same methodology to the binary sample is shown in 
Figure 3b. It is clearly seen that in this case we have a complete sample up to Bt — 13.5. For 
m;i ra > 13.5 we start losing pairs in our sample because it is not possible to select primaries 
with companions fainter than the limiting magnitude of the catalogue. In Figure 3b the solid 
triangles refer to the results of the (V/Vmax) test taking all the binaries for which redshift is 
available for at least one component and we then assumed that both are at the same distance. 
We also show in Figure 3b, as a solid square the (V/Vmax) test taking only those pairs for 
which there is redshift available for both components and AV < lOOOkm/s. 

3 Luminosity Function of Binary Galaxies 

The luminosity function (LF) was determined following a similar procedure used by 
Mendes de Oliveira & Hickson (1991) for compact groups. For a magnitude limited sample 
we have to normalize the contribution of each pair by its effective volume which is defined in 
a such way that the binary put at a given maximum distance would still be included in the 
sample. We have used the same completeness function as used by Hickson, Kindl & Auman 
(1989). Therefore, the expressions used to estimate the effective volume and the LF itself 
were those presented in Mendes de Oliveira & Hickson (1991). As discussed in the previous 
section, using the (V/V max ) test we decided to use a limiting magnitude of 14.75. Also we 
restricted our samples to those pairs for which we know the redshift of at least one member 
and then assumed that both galaxies are at the same distance. Only pairs with galactic 
latitude | b |>30° were included in the LF estimation and the galactic extinction correction 
was assumed to be a 0.12 cossec | b | law. The interval of absolute magnitude adopted to 
estimate the LF was 0.5. 

Figure 4 shows the LF of the binary galaxies in our sample. For comparison we also 
plot the LF obtained by Xu & Sulentic (1991, hereafter XS). It is important to remind that 
these authors measured the LF in a different passband and using a different procedure. Apart 
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of an offset it is quite remarkable the agreement, even confirming a small dip for Mb between 
—20 and —18. We fit the LF with the Schechter function using a nonlinear least-square fit 
program (GAUSSFIT) described in Jefferys et al. (1988). The best fit values for the three 
parameters a, <f>*, and M* are -1.47 ±0.16, 4.8 ±2.2 x 10~ 5 , and -21.23 ±0.20 taking all the 
observed points of the LF into account; and —1.36 ± 0.07, 7. Oil. 7 x 10 -5 , and —21.07 i 0.16 
excluding the two brightest points of the LF. We also superposed on this plot the LF for 
field galaxies obtained by Efstathiou, Ellis, & Peterson (1988) with a = — 1.07 i 0.05, <f>* = 
1.95i0.34 X 10 -3 , and M* = —21.18 i 0.10. It is interesting to note that paired galaxies 
represent ~10% of field galaxies, corroborating previous results from XS. However, while XS 
found that this percentage is approximately constant over the entire luminosity range we can 
only confirm such trend for pairs brighter than Mb = — 19. We should emphasize that both 
works suffer from serious incompleteness problems at the faint end of the LF, which in part 
can explain the discrepancy. 

An estimate of the luminosity density of galaxies in binaries can be done using a, <f>* 
and M* obtained by fitting the LF. Using the parameters obtained with all the observed 
points of the LF we estimate C = L* <f>* T(a ± 2) = 3.9 xlO 6 ~Lq Mpc -3 . The luminosity 
density of field galaxies, as estimated by Efstathiou, Ellis, & Peterson (1988), is 0.96 xlO 8 
~Lq Mpc -3 , implying that binary galaxies represent only 4% of the field galaxies. This low 
percentage compared to the value obtained by XS it is probably due to the fact that we are 
taking a different field LF. 

4 Discussion 

We have used the ESO-LV to build a homogeneous sample of binary systems in the 
southern sky. The catalogue consists of 621 double galaxies, subsequently reduced, by means 
of more restrictive criteria, to a sample of 189 isolated pair of galaxies which are presented 
in tabular form with the galaxy relevant parameters. This is the first objective catalog of 
binary systems in the southern sky and its general characteristics are very similar to those 
previously found for CPG and CMG. 
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The mean separation of pairs in our sample is about 6.9 times the mean diameter of 
both components, and 4% of the pairs in our final sample have the ratio of diameters larger 
than 2. This results is quite comparable to what is found in other catalogues (Sulentic 1992). 
Using those pairs with at least one radial velocity we estimate that their mean separation 
is 143 Kpc. Therefore our list contains contains more wider separation pairs than the CPG 
where the mean separation is 83 Kpc, and is more comparable with the CMG where the mean 
separation is about 146 Kpc (Soares, 1989). 

The mean difference in apparent magnitude between components of CPG pairs is 
0.7±0.3 while in our data this quantity is 0.97 ± 0.67. We conclude that algoritm has iden- 
tified pairs with slightly larger difference in magnitude between the primary and secondary 
objects, although the difference is within \-o . 

About 19% of objects in our sample are E or SO, while 81% are spirals. Specifically, our 
sample have 3% of EE pairs, 32% of ES pairs and 65% of SS pairs, close to what is expected 
from a sample where the objects are randomly associated, which gives 4%EE, 30%ES and 
66%SS pairs. Therefore, our data does not support any evidence of overrepresentation of 
pairs with concordant morfological type. 

We are presently conducting a redshift survey of galaxies in our sample in order to 
produce a homogeneous sample of pairs that allow us to study the physics of interactions 
based on a more representative sample than what it is presently available. 

Acknowledgments — We thank Dr. Steve Zepf for useful suggestions. This work was sup- 
ported in part by FAPESP (RRdC) Proc. No. 92/2686-9, and CNPq (DSLS) Proc. No. 
300193/90-4. The list of objects in Table I is available via an anonymous FTP copy at CDS. 
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APPENDIX 

The CMG Method of Finding Bound Pairs of Galaxies 



In this appendix we reproduce the basic of the method devised by T.S. van Albada (see 
Soares 1989, chapter 2, by van Albada and Soares) to identify pairs of galaxies which exhibit 
physical association. 

The probability of finding the first neighbour of a given galaxy, in the plane of the sky, 
for an ensemble of non-physically associated galaxies, is 



P^dB = e-* B n 2ir0nd0 , (Al) 

where n is the surface density of galaxies. From eq. (Al), one can calculate the expected 
mean value of the nearest neighbour projected separation, i.e., (#1) = l/(2y/n). Defining the 
auxiliar variable 



1 = W) (A2) 

and inserting it in eq. (Al) we find the frequency distribution of x: 

_ 1 2 7T 

pi(x)dx = e 47rx —xdx , (A3) 

which is independent of the surface density of galaxies. 

The method of identifying bound pairs consists in comparing the actual distribution of 
x, p (x) } with pi{x). At small x one expects to find an excess of p over p\ corresponding to 
physically associated pairs, while at large x, the distributions should coincide. 

The surface density of galaxies in the neighbourhood of a given galaxy is evaluated from 
the expected mean separation to a neighbour of a higher order. For the tenth neighbour, 
(#io ) = r(10.5)/(9!y / 7m), and one can approximate 
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T(10.5) 1.76 

The normalized projected separation to the nearest neighbour can now be calculated 
from 6 1 and 9\q\ 

a 

x = Q x 2^ = 3.52-^. (A5) 

#10 



13 



References 

Arp, H., & Madore, B., 1985, A Catalog of Southern Peculiar Galaxies and Associations, 
Cambridge University Press, Cambridge 

Balcells, M., Borne, K.D., & Hoessel, J.G. 1989 ApJ, 336, 655 

Barnes, J.E. 1988 ApJ, 331, 699 

Borne, K.D., 1988 ApJ, 330, 38 

Davis, M., & Huchra, J., 1982 ApJ, 254, 437 

Efstathiou, C, Ellis, R.S., & Peterson, B.A., 1988 MNRAS, 232, 431 

Hickson, P., Kindl, E. & Auman, J.R., 1989 ApJS, 70, 68 

Jefferys, W.H., Fitzpatrick, M.J. & McArthur, B.E., 1988 Cel. Mech., 41, 39 

Karachentsev, I.D., 1972 Comm. Spec. Astrophys. Obs., 7, 1 

Karachentsev, I.D., 1985 Sov Astron, 29, 243 

Karachentsev, I.D., 1987 Double Galaxies, Nauka, Moscou 

Keel, W.C., & van Soest, E.T.M., 1992 A&AS, 94, 553 

Lauberts, A., & Valentijn, E.A., 1989, The Surface Photometry Catalogue of the ESO-Uppsala 
Galaxies, European Southern Observatory, Munich 

Madejsky, R., 1992 A&A, 265, 478 

Mamon, G.A., 1993 in Gravitational Dynamics and the N-body Problem, ed. F. Combes & 
E. Athanassoula, in press 

Mendes de Oliveira, CM., & Hickson, P., 1991 ApJ, 380, 30 

Nilson, P., 1973, Uppsala General Catalogue of Galaxies, Uppsala Astron. Obs. Ann. 6 
Oosterloo, T.A., 1988, Ph.D. Thesis, University of Groningen 
Rampazzo, R., & Sulentic, J.W., 1992 A&A, 259, 43 
Schmidt, M., 1968 ApJ, 151, 393 
Schweizer, L.Y., 1987 ApJS, 64, 427 

14 



Soares, D.S.L., 1989, Ph.D. Thesis, University of Groningen 

Sulentic, J., 1992, in Morphological and Physical Classification of Galaxies, ed. G. Busarello, 
M. Capaccioli, G. Longo, Kluwer, Dordrecht 

Toomre, A., & Toomre, J., 1972 ApJ, 178, 623 

Turner, E.L., 1976 ApJ, 208, 304 

van Albada, T.S., & Soares, D.S.L., 1993, in preparation 

van Moorsel, G.A., 1982, Ph.D. Thesis, University of Groningen 

van Moorsel, G.A., 1987 A&A, 176, 13 

Verner, D.A., & Chernin, A.D., 1987 Sov Astron, 31, 127 

Xu, C., & Sulentic, J., 1991 ApJ, 374, 407 

Zhenlong, Z., Jiansheng, C., Xiaoying, T., & Yulin, B., 1989 Publ Beijing Astr Obs, 12, 8 



15 



Figure Captions: 

Figure 1 : We plot the excess of probability association (p) versus the dimensionless distance 
(x), discriminated for different morphological types. The solid line is the fitting function for 
x<0.75. We can observe that the excess of probability is almost independent of the magnitude 
difference between the primary and secondary component. The symbols to represent such 
difference were used as follows: Solid circle - A J = 0; Solid square - A J = 1; Solid triangle - 
AJ = 2; Open circle - AJ = 3; Open square - AJ = 4; and Open triangle - AJ = 5. 

Figure 2 : Separation in arcsec versus the magnitude of the secondary (m sec ). In a complete 
catalog we expect that all pairs should be located below the continous line. For those pairs 
with m sec fainter than 14.5 there is a tendency for higher separations, which is a natural bias 
produced by the incompleteness of our catalog. 

Figure 3 : The < V/V max > tests for the pairs in our sample. Panel (a) refers to a sample of 
field galaxies extracted from the ESO-LV catalog, with filled and open symbols representing 
magnitudes corrected and not corrected for Galactic extinction, respectively. Panel (b) shows 
the same test applied to our sample of binary galaxies. Solid squares represent those pairs for 
which redshfit is available for both components and solid triangles when redshift is available 
for at least one component. The dotted line indicates the value of < V/V max > for a complete 
sample. 

Figure 4 : The optical luminosity function of the binary galaxies in our sample (solid circle). 
Open circle exhibit the luminosity function obtained by Xu & Sulentic (1991). We also plot 
the LF for field galaxies (dashed line) as provided by Efstathiou, Ellis, & Peterson (1988). The 
solid line shows the best fitting Schechter function obtained using all the available points. The 
dotted line is the best fitting line when we remove the two brightest points of the luminosity 
function. 
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